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Abstract Measurements of magnetic fields and electric currents in the pre- 
emptive corona are crucial to study solar eruptive phenomena, like flare and 
coronal mass ejections(CMEs). However, spectro-polarimetric measurements of 
certain photospheric lines permit a determination of the vector magnetic field 
at the photosphere. Thus, substantial collection of magnetograms relate to the 
photospheric surface field only. Therefore, there is considerable interest in ac- 
curate modeling of the solar coronal magnetic field using photospheric vec- 
tor magnetograms as boundary data. This numerical modeling is carried out 
by applying state-of-the-art nonlinear force-free field (NLFFF) reconstruction. 
Cartesian nonlinear force-free field (NLFFF) codes are not well suited for larger 
domains, since the spherical nature of the solar surface cannot be neglected 
when the field of view is large. One of the most significant results of Solar 
Dynamic Observatory (SDO) mission to date has been repeated observations of 
large, almost global scale events in which large scale connection between active 
regions may play fundamental role. Therefore, it appears prudent to implement 
a NLFFF procedure in spherical geometry for use when large scale boundary 
data are available, such as from the Helioseismic and Magnetic Imager (HMI) 
on board SDO. In this work, we model the coronal magnetic field above multiple 
active regions with the help of a potential field and a NLFFF extrapolation codes 
in a full-disk using HMI data as a boundary conditions. We compare projections 
of the resulting magnetic field lines solutions with full-disk coronal images from 
the Atmospheric Imaging Assembly (SDO/AIA) for both models. This study 



^ Department of Physics, Drexel University, Philadelphia, 
PA 19104-2875, U.S.A, email: 
tasfawOeinstein. physics. d rexei.eduj email: 
Kevin.M.Olson@drexel.edu' 

^ Max Planck Institut fiir Sonnensystemforschung, 
Max-Planck Str. 2, D-37191 Katlenburg-Lindau, Germany, 
email: wiegelmann@mps.mpg.de, email: inhester@mps.mpg.de 

4 NASA, Goddard Space Flight Center, Code 674, Greenbelt, 
MD 20771, U.S.A. email: peter.j.macneice@nasa.gov 

5 Nati onal Solar Observa tory, Sunspot, NM 88349, U.S.A. 
email: |apevtsov@nso.edLj] 



SOLA: solar_global.tex; 1 January 2013; 2:59; p. 1 



T. Tadesse ot al. 



has found that the NLFFF model reconstructs the magnetic configuration better 
than the potential field model. We have concluded that much of trans-equatorial 
loops connecting the two solar hemispheres are current-free. 

Keywords: Active Regions, Magnetic Fields; Active Regions, Models; Magnetic 
fields. Corona; Magnetic fields, Photosphere; Magnetic fields, Models 



1. Introduction 



Magnetic fields are believed to play a dominant role for active phenomena carried 
out in the solar corona. The study of solar eruptive phenomena requires that we 
understand how magnetic energy is stored in the pre-eruptive corona. One of the 
key questions is: What is the three dimensional (3D) structure of magnetic fields 
and electric currents in the pre-eruptive corona, and how much free energy is 
stored in the field? ( Schrijver and Title, 2005] [Jiang and Feng, 2012 ). To answer 
this question we must measure the magnetic field in the coronal volume. How- 
ever, spectro-polarinietric measurements of certain photospheric lines permit a 
determination of the vector magnetic field at the photosphere. Thus, substantial 
collection of magnetograms relate to the photospheric surface field only. Al- 
though measurement of magnetic fields in the chromosphere and the corona has 



considerably improved in recent decades (Lin, Penn, and Tomczyk, 2000 Liu and 



Lin, 2008), further developments are needed before accurate data are routinely 



available. The problem of measuring the coronal field and its embedded electrical 
currents thus leads us to use numerical modeling to infer the field strength in 
the higher layers of the solar atmosphere from the measured photospheric field. 
Due to the low value of the plasma /3 (the ratio of gas pressure to magnetic 
pressure), the solar corona is magnetically dominated (Gary, 2001). To describe 
the equilibrium structure of the static coronal magnetic field when non-magnetic 
forces are negligible, the force-free assumption is appropriate: 



(V X B) X B = 



(1) 



V • B 



(2) 



subject to the boundary condition 



B = B^i^g on photosphere 



(3) 



where B is the magnetic field and B^j-^g is measured vector field on the photo- 
sphere. Equation ([I]) states that the Lorentz force vanishes (as a consequence of 
J II B, where J is the electric current density) and Equation ^ describes the 
absence of magnetic monopoles. 

As an alternative to real measurements, force-free extrapolation of photo- 
sphericmagnetic fields is currently being used as a the primary tool for the 



modeling of coronal magnetic fields (Inhester and Wiegelmann, 2006 Valori, 



Kliem, and Keppens, 2005 Wiegelmann, 2004| Wheatland, 2004 Wheatland 



and Regnier, 2009 Tadesse, Wiegelmann, and Inhester, 2009 Wheatland and 
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Leka, 2011 Amari and Aly, 2010 Wiegelmann, T. and Thalmann, J. K. and 



Inhester, B. and Tadesse, T. et al, 2012 Jiang and Feng, 2012 Jiang, Feng, and 



Xiang, 2012 Aschwanden and Malanushenko, 2012 Wiegelmann and Sakurai, 



2012|). 

Cartesian nonlinear force-free field (NLFFF) codes are not well suited for 
larger domains, since the spherical nature of the solar surface cannot be neglected 



when the field of view is large ( Wiegelmann, 2007| Tadesse, Wiegelmann, and In- 


hester, 2009 


Tadesse et al, 2012a 


Quo et al, 2012 


). One of the most significant 



results of Solar Dynamic Observatory (SDO) mission to date has been repeated 
observations of large, almost global scale events in which large scale connection 
between active regions may play fundamental role. Therefore, it appears prudent 
to implement a NLFFF procedure in spherical geometry for use when large scale 
boundary data are available, such as from the Helioseismic and Magnetic Imager 



(HMI) on board SDO. DeRosa, M. L. and Schrijver, C. J. and Barnes, G. et al. 



(2009) has studied that different NLFFF models have markedly different field 
line configurations and provided widely varying estimates of the magnetic free 
energy in the coronal volume. The main reasons for that problem are (1) the 
forces acting on the field within the photosphere, (2) the uncertainties on vector- 
field measurements, particularly on the transverse component, and (3) the large 
domain that needs to be modeled to capture the connections of an active region 



surroundings ( 


Wiegelmann, 2007 


Tadesse et o/., 2012c 


). In this s 



points explicitly into account for modeling coronal magnetic field in a full-disk. 

In this work, we apply our spherical NLFFF procedure to a group of active 
regions observed on November 09 2011 around 17:45UT by SDO/HMI. During 
this observation, there were four active regions (ARs 11338, 11339, 11341 and 
11342) along with other smaller sunspots spreading on the disk. We compare 
the extrapolated (both potential and NLFFF) magnetic loops with extreme 
ultraviolet (EUV) observations by the Atmospheric Imaging Assembly (AIA) 



(Lemen, J. R. and Title, A. M. and Akin, D. J. et al, 2012) on board SDO. 



During comparison, we check whether the NLFFF model reconstructs the mag- 
netic configuration better than the potential field model. This comparison can 
be used to evaluate how well the model field lines approximate the observed 
coronal loops. 



2. Instrumentation and data set 



The Helioseismic and Magnetic Imager (HMI) (Schou, J. and Scherrer, P. H. 



and Bush, R. I. et al, 2012) is part of the Solar Dynamics Observatory (SDO) 
and observes the full Sun at six wavelengths and full Stokes profile in the Fe 
slowromancapi@ 617.3 nm spectral line. It consists of a refracting telescope, 
a polarization selector, an image stabilization system, a narrow band tunable 
filter and two 4096 pixel CCD cameras with mechanical shutters and control 
electronics. 

The transverse components of vector magnetic fields suffer from the so-called 
180° ambiguity. The 180° ambiguity for the HMI data in this study has been 
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resolved by an improved version of the minimum energy method 


Metcalf, 1994 


Metcalf, T. R. and Leka, K. D. and Barnes, G. et al, 2006 


Leka, Barnes, and 


Crouch, 2009 


). As described in Leka, Barnes, and Crouch 


(2009 


, in weak-field 



areas, the minimization may not return a good solution due to large noise. The 
noise level is w lOG and w lOOG for the longitudinal and transverse magnetic 
field, respectively. Therefore, in order to get a spatially smooth solution in weak- 
field areas, we divide the magnetic field into two regions, i.e., strong-field region 
and weak-field region, which is defined to be where the field strength is below 
200 G at the disk center, and 400 G on the limb. The values vary linearly with 
distance from the center to the limb. The ambiguity solution in the strong-field 
area is derived by the annealing in the minimization method and released in 
the data series "hmi.B_720s_el5wl332_cutout". Magnetic field azimuths in the 
weak-field area are finally determined by the potential-field acute-angle method 
after the annealing. This is different from the original minimum energy method, 
which uses a neighboring-pixel acute-angle algorithm to revisit the weak field 
area. Figure |3^) shows the vector magnetic fields observed by HMI on November 
09 2011 around 17:45UT. 



3. Optimization principle and preprocessing of HMI data 



Molodensky ( 1969 1 and Aly ( 1989 ) pointed out that vector magnetic fields 
on a closed surface that fully encloses any force-free domain have to satisfy 
the force-free and torque-free conditions. To serve as suitable lower boundary 
condition for a force-free modeling, on average a net tangential force acting 
on the boundary and shear stresses along axes lying on the boundary have 
to reduce to zero (Molodensky, 1969} Aly, 1989; Wiegelmann, Inhester, and 



Sakurai, 2006 Tadesse, Wiegelmann, and Inhester, 2009). Preprocessing method 
as implemented in IWiegelmann, Inhester, and Sakurai ( 2006 ) has to be used to 



fulfill those conditions. The preprocessing scheme of Tadesse, Wiegelmann, and] 
Inhester ( 2009 ) involves minimizing a two-dimensional functional of quadratic 



form in spherical geometry similar to 

B = argmin(Lp) , 



(4) 



where B is the preprocessed surface magnetic field from the input observed 
field Bofcs. Each of the constraints L„ is weighted by an as yet undetermined 
factor fin- The first term (n = 1) corresponds to the force-balance condition, the 
next (n = 2) to the torque-free condition, and the last term (n — 4) controls 
the smoothing. The explicit form of Li, L2, and L4 can be found in [Tadesse, | 
Wiegelmann, and Inhester (20091. The term {n = 3) controls the difference 



between measured and preprocessed vector fields. After preprocessing the HMI 



(Wheatland, Sturrock, and Roumeliotis, 2000 


Wiegelmann, 2004 


) in spherical 


geometry ( 


Wiegelmann, 2007 


Tadesse, Wiegelmann, and Inhester, 2009 


Tadesse 


et al, 2012c 


)■ 
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4. Potential field model 



The simplest way to model the coronal field is to assume that it is potential, 
i.e. that it carries no electric current. Solutions for this model in plane geometry 



have been obtained by Schmidt ( 1964 ), for the case where the vertical component 



of the field is specified at the photospheric boundary. Potential field model 
has now led to an almost routine type reconstruction, used for observational 



purposes (Sakurai, 1989), but also for building initial conditions for dynamical 



MHD numerical simulations (Amari et al., 1996). This assumption has proven to 



be adequate for many quiescent, old active regions and even for the non-eruptive 



Rosa, 2003). 



global coronal-heliospheric interface (Wang and Sheeley, 1990 Schrijver and De 



Studies of the coupling of the coronal field into the heliosphere suggest that 
the global coronal magnetic field is often largely potential. For the practical 
calculation of the global field, the so-called source-surface model has been in- 



troduced ( Schatten, Wilcox, and Ness, 1969j , in which the influence of the solar 
wind is artificially taken account of by the requirement that the field be radial 
at some exterior spherical (source) surface typically at = 2.5i?o from the 
suns center. The potential-field source surface (PFSS) model, uses this concept 
to extrapolate the line-of-sight surface magnetic field through the corona with 
the boundary assumed to be at the source surface. 

The potential field approximation has often been used to deduce the magnetic 
structure of the solar corona from measurements of the photospheric field. It 
obeys the equation V x B = 0, so that the magnetic field can be expressed 
as the gradient of a scalar potential, i.e. B — — V$. Since the magnetic field 
is divergence-free (V • B = 0), the scalar potential obeys the Laplace equation 
V^<I> = 0. Together with the photospheric boundary condition, which is tradi- 
tionally provided by a map of the radial magnetic field, the Laplace equation 
can be solved in the spherical coordinate system (r, 0, 0), where 9 stands for 
colatitude. Using the radial magnetic field observed at Rq and the source sur- 
face assumption at R^ as the boundary condition, together with the spherical 
harmonic expansion in the domain Rq < r < i?^ , we finally obtain the potential 
field model. 



5. Results 

In this study, we have used potential field model and our spherical NLFFF 
optimization procedure to a group of active regions observed on November 09 
2011 around 17:45UT by SDO/HMI instrument. During this observation there 
were four active regions (ARs 11338, 11339, 11341 and 11342, see Fig. [2^) along 
with other smaller sunspots spreading on the disk. In order to accommodate 
the connectivity between those ARs and their surroundings, and to study trans- 
equatorial loops (loops connecting north and south hemisphere), we adopt a non 
uniform spherical grid r, 9, 4> with Ur = 225, = 375, ncf, = 425 grid points in 
the direction of radius, latitude, and longitude, respectively, with the field of view 
of ['^min = 1^0 : '^niax = 2Rq] x [9^^^ = -50° : 0max = 50°] x [c^^^^ = 90° : 
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'/'max — 270°]. The main purpose of this work is to study the magnetic and 
electric current connectivities of the northern and southern solar hemispheres 
and to compare which of the two models (potential or NLFFF model) is best 
suited to fully describe those connectivities in full-disk environment. 

We use our preprocessing routine in spherical geometry to derive suitable 
boundary conditions for force-free modeling from the measured photospheric 



data (Wiegelmann, Inhester, and Sakurai, 2006 Tadesse, Wiegelmann, and In 



hester, 2009). We used a standard preprocessing parameter set Mi = A*2 = 1, 



yLt3 = 0.001 and /i4 = 0.05, which are similar to the values calculated from 



vector data used in previous studies (Wiegelmann, T. and Thalmann, J. K. and 


Inhester, B. and Tadesse, T. et al, 2012 


Tadesse et al, 2012c 


) for HMI data. 



The preprocessing influences the structure of the magnetic vector data. It does 
not only smooths Bt (transverse field) but also alters its values in order to reduce 
the net force and torque. The change in is more pronounced than the radial 
component B^ (radial field) since Bt is measured with lower accuracy than the 
longitudinal magnetic field. 

For computing the potential field, we use the preprocessed radial component 
Br of the HMI-data using spherical harmonic expansion. This potential field 
solution has been used to initialize our NLFFF code. For nonlinear force-free 
fields we minimize the functional Eq. (??). We implement the new term Lphoto 
in Eq. (??) to work with boundary data of different noise levels and qualities 



(Wiegelmann and Inhester, 2010 Tadesse et al., 2011 ). For those pixels, for which 



Bobs was successfully inverted, we allow deviations between the model field B 
and the input fields observed Bobs surface field using Eq. (??), so that the model 
field can be iterated closer to a force-free solution even if the observations are 
inconsistent. In order to control the speed with which the lower boundary is 
injected during the NLFFF extrapolation, we have used the Langrangian multi- 



plier oiv — 0.001 as suggested by Tadesse et al. (2012c I. For more details of the 



method used in this work we direct the readers to the study by [Tadesse et al. 



(2011) 



We plot the surface potential field and NLFFF solutions in Figure [T] for 
selected region excluding the polar regions. In this figure, one can see that 
there are substantial differences between the surface radial and transverse field 
components of the two models. To quantify the degree of disagreement between 
the respective vector components of potential and NLFFF model solutions on 
the bottom surface, we calculate their pixel-wise correlations. The correlation 



were calculated( Schrijver, C. J. and Derosa, M. L. and Metcalf, T. R. et al., 



2006 Metcalf, T. R. and De Rosa, M. L. and Schrijver, C. J. et al, 20081 from 



where and are the vectors at each grid point i on the bottom surface. 
If the vector fields are identical, then C^ec — 1; if _L , then C„ec = 0. 
Table [l] shows the correlation (C^ec) of the 2D surface magnetic field vectors 
of potential and NLFFF models for the radial and transverse components. The 
vector correlation between Bi of the potential and NLFFF surface vector maps 
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Figure 1. Magnetic vector maps of HMI data on part of the lower boundary. The color coding 
shows Br on the photosphere and the black arrow indicates the transverse components of the 
field. The vertical and horizontal axes show latitude, 9 and longitude, (f> on the photosphere 
respectively. 
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Table 1. The correlations between 
the components of surface fields 
from potential and NLFFF models. 





V 


U 


C'vcc 




(Bjv_lfff)i- 
(Bjvlfff)* 


0.937 
0.813 



is much more less than that of respective components. This is due to the fact 
that the transverse components (B^) of the potential field model are computed 
using spherical harmonics expansion from the radial components of the mag- 
netic field used as a boundary condition for NLFFF model. The average vector 
correlations between the radial and transverse components of the two models 
in Table [T] indicates that the NLFFF solutions are somewhat far from potential 
as they carry electric currents. In order to compare our reconstructions with 
observation, we plot the selected fieldlines of the potential and NLFFF models in 
Figure[2]; and d. We overlay the field lines with AIA 193 A image. For observation 
time in Figure ^ and d, the field lines of the potential and NLFFF models are 
reconstructed from the same footpoints. The potential field lines in Figure ^c) 
have an obvious deviation from the observed EUV loops, since the projection 
of the field lines at the bottom of AR 11339 leans toward the east, while the 
loops toward the west. The loops connecting positive and negative polarities of 
AR 11341 are best overlaid by NLFFF lines than potential ones. The spatial 
correspondence between the overall shape of the NLFFF field lines and the 
EUV loops is much improved as shown in Figure [2|d). Therefore, the qualitative 
comparison between the model magnetic field lines and the observed EUV loops 
indicates that the NLFFF model provides a more consistent field for full-disk 
magnetic field reconstruction. In the absence of a more reliable quantitative 
comparison, it remains the best option even if the qualitative nature is not 
ideal. 

In addition to the above comparison to quantify the degree of disagreement 
between vector field solutions of the two models in the computational volume 
that are specified on identical sets of grid points, we use five metrics that compare 
either local characteristics (e.g., vector magnitudes and directions at each point) 
or the global energy content in addition to the force and divergence integrals as 



defined in [Schrijver, C. J. and Derosa, M. L. and Metcalf, T. R. et aZ.| ( |2006 ). 
The vector correlation (C^ec) metric of Equation (|5| is also used analogous to 
the standard correlation coefficient for scalar functions. The second metric, Ccs 
is based on the Cauchy-Schwarz inequality(|a • b < |a||b| for any vector a and 
b) 

where B = Bp^^, b ~ -^NLFFF' ^^'^ ^ ^^"^ number of vectors in the field. 
This metric is mostly a measure of the angular differences of the vector fields: 
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(a) SDO/HAli luj.niLotogram (!>) 17lA AIA iiiiafje 















1 






1 



Figure 2. a) Full-disk SDO/HMI magnctogram, the rectangles figure outline the main active 
regions during this observation and the yellow line indicates the solar equator (Q = 0°), b) 
Full-disk AIA 171 A image. Both data sets were obtained on November 09, 2011 around 
17:45UT, c) Field lines of the potential field model at 17:45 UT overlaid on AIA 193A image 
and d) Field lines of the NLFFF model at 17:45 UT overlaid on AIA 193A image. Green and 
red lines represent open and closed magnetic field lines. 



Ccs — 1 when B and b are parallel and Ccs = — 1 if they are anti-parallel; 
C'cs = if Bi _L bi at each point. Next, we introduce two measures for the 
vector errors, one normalized to the average vector norm, one averaging over 
relative differences. The normalized vector error Ej^ is defined as 

£;N-E|bi-Bi|/E|Bi|, (7) 

i i 

The mean vector error Em is defined as 



1 ^ |bi-Bi| 
nI^ ir,i • 
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Table 2. Comparision of figure of merits between potential field and NLFFF models. We 
have used phorical grids of 225 X 375 X 425. 





Model 






C*vcc 


Ccs 


I- En 


I- Em 


e 


Time 


Potential 
NLFFF 


0.000 
0.591 


0.002 
0.997 


1 

0.878 


1 

0.833 




0.735 




0.876 


1 

1.220 


30min 
4h:39min 



Unlike the first two metrics, perfect disagreement of the two vector fields results 
in 1 — E}^[ = 1 — En = 0. As we are also interested in determining how 
well the models estimate the energy contained in the field, we use the total 
magnetic energy in NLFFF model field normalized to the total magnetic energy 
in potential field as a global measure of the quality of the fit: 



E.|B,| 



2 ■ 



e = 1, if there is no difference between the potential field and the nonlinear 
force-free model solutions. The degree of convergence towards a force-free and 
divergence-free model solution can be quantified by the integral measures of 
the Lorentz force and divergence terms in the minimization functional in Equa- 
tion (??), computed over the entire model volume Y. Lj and Ld of Equation 
(??) measure how well the force-free and divergence-free conditions are fulfilled, 
respectively. In Table [2j we provide some quantitative measures to rate the qual- 
ity (figure of merit) of our reconstruction. Column 1 names the corresponding 
models. Columns 2 — 3 show how well the force and solenoidal condition are 
fulfilled for both models. The next five columns of Table |2] contain different 
measures which compare our NLFFF solution with potential field. Those figures 
of merit indicate that there is a clear difference between potential field and 
NLFFF models solutions for full-disk. The last column shows the computing 
time on 1 processor. The figures of merit show that the potential field is far 
away from the true solutions and contains only 81.97% of the total magnetic 
energy. 

We estimate the free magnetic energy to be the difference between the ex- 
trapolated NLFFF and the potential field with the same normal boundary 



conditions in the photosphere(Regnier and Priest, 2007 Thalmann, Wiegelmann, 



and Raouafi, 2008). We therefore estimate the upper limit to the free magnetic 



energy associated with coronal currents of the form 

^frcc ^^J^ {bIi f f - Blot) r''sin9drdedct>, (10) 

Our result for the estimation of free-magnetic energy in Table |3] shows that the 
NLFFF model has 18.83% more energy than the corresponding potential field 
model. Therefore, whenever we try to analysis magnetic field configuration and 
its magnetic field contents in full-disk, it is inaccurate to use potential field model 
as it undermines the real features of the field in the corona. 
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Table 3. The magnetic energy associated with 
extrapolated potential and NLFFF field configu- 
rations from full disk SDO/HMI data. 





Model 


^totalClO^'e'-g) 


%gg{1033erg) 


Potential 
NLFFF 


7.306 

8.913 




1.607 



In our previous work (Tadesse et al., 2012a), we have studied the connectivity 



between three neighbouring active regions. In this work, we study the magnetic 
and electric connectivities between active regions in northern and southern solar 
hemispheres using full-disk HMI data. During this observation, there were three 
active regions in the northern hemisphere and one active region with surrounding 
sunspots in the southern hemisphere. As a result the total unsigned magnetic 
field flux in the northern hemisphere is larger than that of the southern one. In 
order to quantify these connectivities, we have calculated the magnetic flux and 
the electric currents shared between those active regions. For the magnetic flux, 
e.g., we use 

$„/3 = 51 1^* • f\Rlsin{e,)M,^4>, (11) 

i 

where the summation is over all pixels of AR^ from which the field line ends 
in AR^ or ? G ARq|| conjugate footpoint(j) G AR^. The indices a and (3 de- 
note the active region. For the electric current we replace the magnetic field, 
Bi • r, by the vertical current density • f in Equation (111. Table |4] shows 



the percentage of the total magnetic flux and electric current shared between 
the ARs 11339 and 11338 and between AR 11341 and sunspots in the southern 
hemisphere (see Fig. [S]). We have calculated total an unsigned flux for each 
active regions 11339 and 11341 (both being in northern hemisphere) and the 
flux due to those field lines ending in the AR 11338 and magnetic patches in 
the southern hemisphere, respectively. The first column of Table [3] shows that 
7.31%(6.17%) of positive/negative polarity of AR 11339 (11341) in the northern 
hemisphere is connected to positive/negative polarity of AR 11338 (magnetic 
patches) in southern hemisphere for potential field configuration. The percentage 
share between those ARs are 7.57%(6.91%) for NLFFF configuration. In Table|4) 
the percentages in the bracket show the share in electric current between the 
corresponding active regions and small magnetic patches. Figure |3] display the 
trans-equatorial loops connecting those active regions. In the figure, we rotate 
those active regions to the limb to show loops connecting the two hemispheres. 
The surface contour plots in Figures |4] and [5] show the magnetic flux and electric 
current density flux of field lines crossing equatorial plane, respectively. From 
those figures we can see that NLFFF model has more magnetic and electric 
current density fluxes compared to the potential field. Figure [5] shows that there 
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are more electric currents crossing the plane carried by trans-equatorial loops 
connecting AR 11341 in northern hemisphere to sunspots in the south than 
those currents connecting ARs 11339 and 11338. This indicate that much of the 
trans-equatorial loops connecting ARs 11339 and 11338 are potential. In general 
the two hemispheres are more magnetically connected than electric current. 
However, one has to use NLFFF model over potential model to study magnetic 
field structure in any phenomena involving the field. 




(c) Potentuil fifid (d) NLFF iicid 



Figure 3. a) Selected potential field lines connecting ARs 11338 and 11339 rotated to limb, 
b) The corresponding field lines for NLFFF in a), c) Selected potential field lines connecting 
AR 11341 and Sun spots in southern hemisphere rotated to limb d) The corresponding field 
lines for NLFFF in c). 



6. Conclusion and outlook 

Both potential and nonlinear force-free field (NLFFF) codes in Cartesian ge- 
ometry are not well suited for larger domains, since the spherical nature of the 
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(b) M.FK fioid 



j/(gcu3s) 



-!UU -dU -■iU -4 LI -liiO 20 40 bO bU- 100 

Figure 4. Surface contour plot of perpendicular component of magnetic fields crossing a plane 
of the solar equator (6 = 0°). The color coding shows _Bx crossing the 5 = 0° equatorial plane. 
The vertical and horizontal axes show height, r(in solar radius) and longitude, (/)(in degree) on 
the plane. 
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Figure 5. Surface contour plot of perpendicular component of electric current density crossing 
a plane of the solar equator (6 = 0°). The color coding shows J± crossing the 9 = 0° equatorial 
plane. The vertical and horizontal axes show height, r(in solar radius) and longitude, 4){m 
degree) on the plane. 
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Table 4. The percentage of the total magnetic flux and electric current 
density shared between the ARs 11339 and 11338 and between AR 11341 
and sunspots in the southern hemisphere ( see Fig. |3]l. The percentage in 
brackets denote that of electric currents. 





Model 


AR 11339 1 — > AR 11338 


AR 11341 1 — 5> Patches in South 


Potential 
NLFFF 


7.31 (0.0) 
7.57(0.13) 


6.17(0.01) 
6.91 (3.35) 



solar surface cannot be neglected when the field of view is large. One of the 
most significant results of Solar Dynamic Obseratory (SDO) mission to date has 
been repeated observations of large, almost global scale events in which large 
scale connection between active regions may play fundamental role. Therefore, 
it appears prudent to implement a NLFFF procedure in spherical geometry for 
use when large scale boundary data are available, such as from the Helioseismic 
and Magnetic Imager (HMI) on board SDO. 

In this study, we have investigated the coronal magnetic field associated with 
full solar disk on November 09 2011 by analysing SDO/HMI data using potential 
and NLFFF models. During this particular observation, there were three active 
regions in the northern hemisphere and one active region surrounded by sunspots 
in the south. We have used our spherical NLFFF and potential codes to compute 
the magnetic field solutions in full-disk. For computing the potential field, we use 
the preprocessed radial component of the HMI-data using spherical harmonic 
expansion. We implement our NLFFF code initialized by the this potential field 
solution (except the lower observed bottom boundary) during relaxation towards 
force-freeness state in the computational volume. 

We have compared the magnetic field solutions from both models. The quali- 
tative comparison between the model magnetic field lines and the observed EUV 
loops indicates that the NLFFF model provides a more consistent field for full- 
disk magnetic field reconstruction. In addition to this, the figures of merits have 
been used to quantify the disagreements between the two models field lines. 
The figures of merit indicate that the magnetic field lines structures are not 
potential in full-disk field-of- views. However, much of the trans-equatorial field 
lines connecting the active regions in the northern and southern hemispheres are 
potential. This indicates that the two solar hemispheres are more magnetically 
connected than electric current. The magnetic field lines obtained from nonlinear 
force-free extrapolation bear larger total magnetic energy that of the potential 
field model. Therefore, one has to use NLFFF model over potential model to 
study magnetic field structure in any phenomena involving the field in full-disk. 
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